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Abstract 

We investigate theoretically a fiber laser passively mode-locked with non- 
linear polarization rotation. A unidirectional ring cavity is considered with 
a polarizer placed between two sets of a halfwave plate and a quarterwave 
plate. A master equation is derived and the stability of the continuous and 
mode-locked solutions is studied. In particular, the effect of the orientation 
of the four phase plates and of the polarizer on the mode-locking regime is 
investigated. 



1 Introduction 



Passively mode-locked fiber lasers are of great importance for various applica- 
tions involving optical telecommunications. Different experimental methods 
have been used to achieve mode- locking operation PQ- JI]. In this paper we 
are interested in mode-locking through nonlinear polarization rotation. This 
technique has been successfully used to obtain short pulse generation in dif- 
ferent rare-earth doped fiber lasers and is self-starting. The 
laser configuration is a unidirectional fiber ring cavity containing a polarizer 
placed between two polarization controllers. The polarization state evolves 
nonlinearly in the fiber as a result of the optical Kerr effect. If the polar- 
ization controllers are suitably oriented, the polarizer lets pass the central 
intense part of a pulse while it blocks the low intensity wings. 

Different theoretical approaches have been developed to describe the 
mode-locking properties of such laser. Haus et al. [TJ |2| have developed a 
model based on the addition of the different effects assuming that all effects 
are small over one round-trip of the cavity. Analytical studies of Akhmediev 
et al. ^21 E] are based on a normalized complex cubic Ginzburg-Landau 
(CGL) equation and give the stability conditions of the mode-locked solu- 
tions. On the other hand, many numerical simulations have been done to 
complete analytic approaches ^HJ-jSOl- We have recently investigated ex- 
perimentally and theoretically the mode-locking properties of an Yb-doped 
double clad fiber laser passively mode-locked through nonlinear polarization 
rotation ^21 12] • The optical configuration was a unidirectional ring cavity 
containing an optical isolator placed between two halfwave plates. Only two 
phase plates were considered for simplicity. The theoretical model reduces 
to a complex cubic Ginzburg-Landau equation whose coefficients explicitly 
depend on the orientation of the phase plates. The model allowed the de- 
scription of both the self-starting mode-locking operation and the operating 
regimes as a function of the orientation of the halfwave plates. The model 
was then adapted to the anomalous dispersion case [22] and to the stretched- 
pulse operation [23] • Although our simplified model is in good agreement 
with the experimental results, a typical experiment includes two polariza- 
tion controllers instead of two halfwave plates. Indeed, mode-locking is more 
easily obtained in the former case because there is more degrees of freedom. 
The aim of this paper is to provide a general model taking into account a 
polarizer and two sets of a halfwave plate and a quarterwave plate. The pa- 
per is organized as follows. In section 2 we derive a propagation equation for 
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a unidirectional ring cavity containing a polarizer placed between two sets 
of a halfwave and a quarterwave plates. The resulting equation, valid for 
a large number of round trips, is of the CGL type and explicitly takes into 
account the orientation of the phase plates and the polarizer. Constant and 
mode-locked solutions are considered in section 3. The last section is devoted 
to a discussion of the influence of the phase plates and of the polarizer on 
the stability of both the mode-locking and the continuous wave regimes of 
the laser. 

A/2 A/4 Polarizing A/4 A/2 
n n isolator n n 




Figure 1: Schematic representation of the fiber laser passively mode-locked through 
nonlinear polarization rotation. 



2 The master equation 

In this section we derive a master equation for the laser shown in figure ^ 
The ytterbium-doped fiber has gain, birefringence, group velocity dispersion 
(GVD) and optical Kerr nonlinearity. The cavity contains a polarizing iso- 
lator placed between two polarization controllers. 

2.1 Propagation along the ytterbium-doped fiber 

In the framework of the eigenaxis of the fiber moving at the group veloc- 
ity, the propagation equations for the two polarization components of the 
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amplitude of the electric field are |T2J 12H HHj 
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where g in m _1 is the linear gain, c<j 9 = 10 13 s _1 the spectral gain bandwidth, 
A = 2/3, and 5 = 1/3. K is the birefringent parameter and 7 the nonlinear 
coefficient. 

Following our analysis of reference jT^j, we assume that the effects of the 
GVD /3 2 , the nonlinear effect 7, and the gain filtering p = g/uJ 2 are small 
over one round-trip of the cavity. A perturbative approach can be used. We 
introduce a small parameter e and replace the quantities fa, 7 and p by efa, 
£7 and ep. Let (u(0),v(0)) the electric field components at the entrance of 
the ytterbium-doped fiber, and (u(L),v(L)) the components at the exit of 
the fiber of length L. A first order perturbative calculation leads to ^21 122j 
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2.2 Modelling the phase plates and the polarizer 

The Jones matrix formalism is well adapted to the treatment of a combination 
of phase plates and polarizer. It will be used in this section. Without loss of 
generality, we assume that the eigenaxis at both ends of the fiber are aligned 
and parallel to the x and j/-axes of the laboratory frame. Let et\ (resp. a^) 
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the angle between the eigenaxis of the halfwave plate and the x-axis before 
(resp. after) the polarizer. Let a 2 (resp. 0:3) the angle between the eigenaxis 
of the quarterwave plate and the x-axis before (resp. after) the polarizer. 
Let 9 the angle between the passing axis of the polarizer and the x-axis. 

In the framework of their eigenaxis, the Jones matrices of the quarterwave 
and halfwave plates are respectively 

-*-2(ViJ«). ^ 
Mo!)- <«» 

Let M 3 (resp. M 4 ) be the Jones matrix of the quarterwave plate (resp. 
halfwave plate) after the isolator in the (Ox, Oy) frame: 

M 3 = R(a 3 )M,R(-a 3 ), (7) 

4 

M 4 = R(a 4 )MxR(-a 4 ), (8) 



where 



, cos a — sin a . 
R(a) = . (9) 
1 sina cos a 



is the rotation matrix of angle a. 

Light exiting the polarizer passes through a set of a quarterwave and a 
halfwave plates. Therefore the electric field at the entrance of the fiber after 
the n th round trip is 

ilSjHMr)' (io) 

where u' n and v' n are the electric field components just after the polarizer. 

Let M be the Jones matrix of the polarizer and Mi (resp. M 2 ) the Jones 
matrix of the halfwave plate (resp. quarterwave) before the polarizer. In the 
(Ox, Oy) frame, the matrices write as 

M = i?w (o (ii) 

where f3 = 95% is the transmission coefficient of the polarizer, and 

Mi = R(a l )MxR{-a l ), M 2 = R{a 2 )MxR{-a 2 ). (12) 
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The field after the polarizer can be written as 



(13) 

where f n+ i is the electric field amplitude after the polarizer at the (n + l) th 
round trip. 

We now replace the matrices M, M±, and M2 by expressions and 
(H2J) respectively. We further take for (u n (L) , v n (L)) the expressions given in 
(|3l4j) . and (w„(0), f„(0)) is replaced by equation (JTUJ). Finally, we take into 
account equations (J7|) and (jSJ), and get a relation between f n+1 and / n : 

/n+i = /3e ffi |g/ n + £ ^ p -^LQ^ + l Pf n \f n \ 2 } + 0( £ 2 ), (14) 

where the coefficient P and Q are given in the appendix. The important fact 
in our analysis is that coefficients P and Q explicitly depend on the angles 
cti, 0:2, 03, ct4, and 9. As we will see in the next section, the model will 
allow to investigate the operating regime of the laser as a function of the 
orientation of the phase plates and of the polarizer. 

A stationary state is reached when |/ n +i| = |/n|- This occurs when the 
gain attains its threshold value g = go + £gi + O (e 2 ). g\ is referred to as the 
excess of linear gain below. The dominant part of f n+ i is obtained at order 



f n+ i = l3e 9oL Qf n + 0{e). (15) 

As a consequence of the stationarity, the modulus of Pe 9oL Q is unity. We 
thus obtain the expression of go, as 

So= ^ln03 2 |Q| 2 ) 

2 _ L (16) 

= 2I ln O^il 2 + <? iKL m + e- 2lKL M* 2 + |0 2 | 2 ]) . 

By performing a Taylor expansion of e £9lL , and replacing (3e 9oL Q by e 1 ^, 
equation ((Hj) becomes 

f n+1 = e** (1 + e 9l L) f n + e(p~^ Le**^ + ^Pf n \f n \ 2 + O (e 2 ) . 

(17) 
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It is more convenient to describe the evolution of the field amplitude f n by a 
continuous equation. The discrete sequence f n is interpolated by a continuous 
function and, for a large number of round trips noc 1/e, a fast rotating phase 
factor is set apart (T21 122] , which yields the equation 

i-QC = igiF +[Y + tp )m^ + {Vr + lVi) F |F|2 ' (18) 

where 

F (C = enL) = f n e~ 1 ^ + O (e) , (19) 
and T> r and Z>j are the real and imaginary parts of the quantity T> given by 

v = wv < 2 °) 

They correspond respectively to the effective self-phase modulation and to 
the effective nonlinear gain or absorption. T> r is always negative while the 
sign of T>i depends on a%, a 2 , «3, «4, and 6. Equation (JTSj) is of cubic complex 
Ginzburg-Landau type (CGL). 

3 Solution of the CGL equation 

This section is devoted to the study of two particular solutions of equa- 
tion We first consider the constant solution corresponding to a contin- 
uous wave (CW) operating regime of the laser. Localized solutions are then 
considered and are related to the mode-locking regime of the laser. In both 
cases, the stability criterium of the solution is determined. 

3.1 Constant amplitude solution 

A constant amplitude solution of CGL is 

F = A e m- nt \ (21) 

where 

Q 2 = - (V,\A\ 2 + 9l ), k = ^ (V,\A\ 2 + 9l ) - V r \A\ 2 . (22) 
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Solution (|2T|) is time independent if fl = 0. Under this condition, the expres- 
sions of A and k are 



This solution exists only if T>ig\ is negative. On the other hand, it has been 
demonstrated that the modulational instability occurs when the excess of 
linear gain g\ is negative and the effective nonlinear gain T>i is positive [T2*] . 
Therefore the constant amplitude solution is stable when the excess of linear 
gain is positive and the effective nonlinear gain Pj is negative. 

3.2 Localized solution 

Equation (j!8|) admits the following localized solution: 

F = a{t) 1+id e~^, (24) 

where 



_ -3 [(3 2 V r + 2pVi] + ^9 [2pV t + (3 2 V r } 2 + 8 - 2pV r ] 2 

d= 2 [p 2 V t - 2 P V T ] ' (25) 

-gi [Apd + p 2 d 2 - (3 2 \ 

U ~ 2{pdi-p-(3 2 d} ■ (26) 

The parameter d represents the chirp. The amplitude a(t) writes as 

a(t) = MNsech(Mt), (27) 

where 

M =y^&- (28 » 



The pulses exist if both M and are real. Stability of the localized solution 
results from an equilibrium between the excess of linear gain, the quantity 
/3 2 V r , and the effective nonlinear gain. Indeed, in the defocusing case where 
(3 2 V r < 0, the pulse is potentially stable if the excess of linear gain g 1 is 
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negative and the effective nonlinear gain T>i is positive. This criterium can 
be written in the mathematical form ^2] 

(pd 2 - p - fad) < 0. (30) 

When the effective nonlinear gain is negative, the stability of the pulses is 
not known at this time. Note that higher order terms or gain saturation can 
definitely stabilize the short pulse solution of equation (flHjl. 

4 Influence of the orientations of the phase 
plates and of the polarizer 

In the previous section we have derived a master equation for a laser passively 
mode-locked by nonlinear polarization rotation. The coefficients of the equa- 
tion depend on the orientation angles of the phase plates a%, a 2 , 03, «4, and 
of the polarizer 9. As a consequence, the stability of both the continuous and 
the mode-locked solutions also depends on these angles. Because of the large 
number of degrees of freedom, we cannot perform a systematic study of the 
stability of the solutions as a function of the five angles. In the following we 
have generally fixed three angles and varied the two remaining ones. In these 
conditions it is convenient to summarize the results in a two dimensional 
stability diagram which gives for any couple of varying angles the regions 
of stability of both the continuous and the mode-locked solutions. We have 
first considered (9, a 2 , ^3) = ($> 0°, 0°) where 9 takes the following values: 0°, 
30°, 45°, 60°, 90°, 120°, 135°, 150°, and 180°. We have plotted the stability 
diagram in the plane (oti, 04) for each value of 9. The same studies have been 
done for (9, a 2 , a 3 ) = (0°, a 2 , 0°), (0°, 0°, a 3 ), (30°, 30°, 30°), (45°, 120°, 150°), 
and (60°, 30°, 135°). In the two first cases, a 2 and a% take the same values 
as attributed to 9. For the numerical computations, we have used the same 
parameters as in ref. ^2] : K = 1.5 m _1 , (3 2 — 0.026 ps 2 m _1 , L = 9 m and 
7 = 3-10- 3 W _1 m- 1 . 

A great dependance of the stability domains versus a±, a 2 , 03, «4, and 
9 have been observed. This can be physically expected because a change 
in the orientation of one element leads to a relative variation of the losses 
undergo by the wings and the center of the pulse. It is then possible either 
to favor the center of the pulse which travels the polarizer with a minimum 
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losses, leading to efficient mode-locking regime, or to favor the opposite case 
resulting in the instability of the mode-locking regime. These results are 
illustrated in figures El El IH and El They give the stability domains 

90° 
a 4 (deg.) 



45° 



0° 

0° 45° 90° 

<*i(deg.) 

Figure 2: Stability diagram of the CW and the mode- locked solutions in the plane 
(0:1,0:4) for (^,02,0:3) = (0°,0°,0°). The white region corresponds to stable CW opera- 
tion and unstable mode-locking, the light gray corresponds to unstable CW and unstable 
mode-locking and the dark gray region corresponds to stable mode-locking operation and 
unstable CW. 

of the CW and mode-locking regimes depending on the orientation angles 
(oti, 0:4) of the halfwave plates, for the following orientations of the polarizer 
and quarterwave plates: (0,a 2 ,a 3 ) = (0°,0°,0°), (0°,0°,30°), (0°,45°,0°), 
and (0°,0°,45°), respectively. The representations have been limited to 0° ^ 
ax, a 4 ^ 90° because of the periodicity. Figure El is the same that the one in 
reference ^2] where only two halfwave plates were considered. This is correct 
because the polarizer is aligned with the eigenaxis of the two quarterwave 
plates. Thus this result validates the general model including four phase 
plates. A large part of the computed cartographies are relatively close to 
figure El but another typical shape is shown on figure El Figures HI and El 
show that the operating regime can be independent of the orientation of one 
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0° 45° 90° 

tfi(deg.) 

Figure 3: Stability diagram of the CW and the mode- locked solutions in the plane 
(0:1,0:4) for (^,02,0:3) = (0°,0°,30°). The white region corresponds to stable CW oper- 
ation and unstable mode-locking, the light gray corresponds to unstable CW and mode- 
locking and the dark gray region corresponds to stable mode-locking operation and unsta- 
ble CW. 



of the halfwave plates. We can note on figure 03 that the orientation of the 
last half wave plate (04) does not modify the stability of the solutions in this 
case. This is not surprising since for 8 = 0° and 03 = 45°, the polarization 
that enters this last half wave plate is circular. Whatever the orientation of 
this plate, the polarization entering the fiber is thus circular, and the global 
behavior does not depend on a 4 . This further allows to give a physical 
interpretation to the absence of any mode-locking domain in this case. We 
can see from relations (3) and (4) that in absence of birefringence (K = 0), 
if a circular polarization enters the fiber, a circular polarization exits the 
fiber. Actually, nonlinear polarization rotation does not occur. We can thus 
assume that this is the reason why no mode-locking regime is predicted here. 

We have then explored the dependency of the operating regimes of the 
laser with respect to the orientation angles (a 2 , a 3 ) of the quarterwave plates. 
The periodicity versus ct 2 and a 3 is 180°. Figures El and |H] give typ- 
ical examples of cartographies. They are obtained for the orientations 
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0° > ■ ^ 

0° 45° 90° 

<*i(deg.) 

Figure 4: Stability diagram of the CW and the modedocked solutions in the plane 
(ai, 014) for (0, Q!2, 0:3) = (0°, 45°, 0°). The colors have the same meaning as in figure|21 



(9, «i, on) = (0°, 0°, 0°), (0°, 30°, 0°), and (60°, 30°, 135°) of the polarizer and 
halfwave plates. We can note the large regions of instability and also the in- 
creased number of mode-locking regions compared to the reference results of 
figure El especially on figure |H1 It is interesting to point out the existence of 
four horizontal axes that separate abruptly the different domains and where 
no mode-locking is observed. They locate at values of 03 about integer mul- 
tiples of 45°, on figures EK and around 15°, 45°, 105° and 135° on figure |H1 
In the latter case, 9 = 60°, while it is zero in the former. We can thus deduce 
that for 0:3 = 60° ± 45°, polarization exiting the plate n°3 is circular, which 
is not modified by the last plate n°4 (A/2). As previously, we can assume 
that nonlinear polarization rotation does not occur such that mode-locking 
is not observed. These cases correspond indeed to the horizontal axes where 
a 3 is around 45° or 135° on figures EHZl 15° or 105° on figure |H1 In addition, 
these axes appear as boundaries: when a 3 passes through these axes, the 
ratio between the x-polarized and the y-polarized components entering the 
fiber passes unity, "inverting" the effect of nonlinear polarization rotation 
and thus on mode-locking or CW operation. We have checked with other 



11 
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a 4 (deg.) 
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90° 



Figure 5: Stability diagram of the CW and the mode- locked solutions in the plane 
(ai, at) for (9, a>2, as) = (0°, 0°, 45°). The colors have the same meaning as in figure^! 



values of 6 the existence of similar horizontal axes at = 6 ± 45° , that 
separate abruptly mode-locking and CW domains and where mode-locking 
does not occur in general. Other axes, around = 0° and 90° on figures H30 
or 45° and 135° on figure |H1 can be interpreted with similar arguments. The 
eigenaxes of this plate are then parallel to those of wave plate n°4 (0:4 = 0° in 
the former case, 135° in the latter). Then the polarization entering the fiber 
is in general elliptical, but with its high-axis oriented at 45° from the x-axis 
and y-axis of the fiber. The maximum of x and y amplitudes in the fiber are 
thus identical and we can assume that nonlinear polarization rotation is not 
efficient. To confirm this assumption, we have plotted another cartography 
in the (^2,03) plane with the same parameters: 6 = 30° and a.\ = 30°, but 
with «4 = 120° (not drawn here). In this case, two horizontal axes without 
any mode-locking are located at a.3 = 15° and 105° instead of 45° and 135°. 
These axes correspond to orientations such that the polarization entering the 
fiber is elliptical with its high-axis oriented at 45° of the x-axis and y-axis of 
the fiber. This is thus similar to previous cases with = 0° or 135° and we 
can understand that no ML occurs for these two horizontal axes. Note that in 
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90° 
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180° 



Figure 6: Stability diagram of the CW and the mode- locked solutions in the plane 
(a2, as) for (9, oc\, 014,) — (0°, 0°, 0°). The colors have the same meaning as in figure|21 



this case, two other axes are observed for near 75° and 165°. Polarization 
exiting the plate n°3 is then circular, which is not modified by the last plate 
n°4. Nonlinear polarization rotation is then very difficult to be obtained, as 
already mentioned. 

We have seen that it is possible to give some physical interpretations 
concerning the influence of parameters a 3 and a 4 , located just before the 
fiber. Polarization states can then be well understood since these elements are 
located just after the polarizer. In contrast, it is very difficult to interpret the 
influence of parameters a\ and «2 located at the exit of the fiber. Influence 
of these parameters depends indeed strongly on polarization effects induced 
in the fiber, which are not directly accessible. Experimentally the role of 
phase plates n°l and n°2 is essential because they allow the adjustment 
of the polarization state of the incident electric field at the entrance of the 
polarizer, in such a way that the central part of the pulse is transmitted while 
the wings are blocked. However no quantitative description of the influence 
of the orientation of phase plates n°l and n°2 has been found, due to the 
high complexity of the nonlinear dynamics. We but point out their key role. 

Let us now consider the influence of the orientation of the polarizer 9 on 
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0° 90° 180° 
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Figure 7: Stability diagram of the CW and the mode-locked solutions in the plane 
(tt2, Qfa) for (8, ot\, a±) = (0°, 30°, 0°). The colors have the same meaning as in figure|21 



the operating regimes of the laser for fixed orientations of the phase plates. 
Some diagrams are represented in figureElfor a 2 , a^, a^) = (0°, 0°, 0°, 0°) 
(a), (30°, 45°, 120°, 150°) (b), and (30°, 0°, 0°, 30°) (c). We can note on these 
figures and also on many diagrams not reported here that for any values 
of the orientations of the phase plates, mode-locking can be achieved by a 
rotation of the polarizer. 

In summary, although some behaviors can be well interpreted, it is very 
difficult to deduce general trends for the mode-locking properties of the laser 
essentially because of the large number of variable parameters. However, the 
model is a very powerful tool to predict the behavior of the laser. 

5 Conclusion 

In conclusion we have developed a general model for a fiber laser passively 
mode-locked by nonlinear polarization rotation. A unidirectional ring cavity 
containing a polarizer placed between two sets of a halfwave and a quarter- 
wave plates each has been considered. Starting from two coupled nonlinear 
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Figure 8: Stability diagram of the CW and the mode- locked solutions in the plane 
(o!2j ots) for (9, a>i, a±) = (60°, 30°, 135°). The colors have the same meaning as in figure[21 
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Figure 9: Stability of the CW and the mode-locked solutions versus 9 for 
(ai,a 2 ,a 3 ,a 4 ) = (0°, 0°, 0°, 0°) (a), (30°, 45°, 120°, 150°) (b), (30°, 0°, 0°, 30°) (c). The 
colors have the same meaning as in figure 



propagation equations for the electric field components we have derived a 
unique equation for the field amplitude, which is a complex cubic Ginzburg 
Landau equation. The coefficients of the equation depend explicitly on the 
orientation angles of the polarizer and of the phase plates. We have thus 
investigated the stability of both the constant amplitude and the short-pulse 
solutions as a function of the angles. Solutions have been found analytically. 
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Although it is difficult to give some general trends, the model has the ad- 
vantage to describe a real experiment. Indeed, it includes the linear and 
nonlinear characteristics of the doped fiber, two polarization controllers and 
a polarizer. 



Appendix 

We give hereafter the coefficients of the master equation: 

Q = e - lKL ^ + e iKL <f) 2 , 

01 = (Xi cos Q + X2 sin 9) (x 3 cos 9 + % 4 sin 9) , 

02 = (X3 sin - xl cos e ) (x*i sin e ~ X*2 cos O) , 
-V2 



Xi 

X2 
X3 
X4 



[(i + cos(2« 3 )) cos(2a 4 ) + sin(2a 3 ) sin(2a;4)] , 
[(i — cos(2a: 3 )) sin(2o;4) + sin(2a 3 ) cos(2o;4)] , 



2 
V2 
2 

-V2 



[(i + cos(2a 2 )) cos(2«i) + sin(2«i) sin(2a 2 )] , 
[(i — cos(2a2)) sin(2«i) + cos(2o;i) sin(2a2)] , 



(31) 

(32) 
(33) 

(34) 
(35) 
(36) 
(37) 



and 
P = e 



with 



-iKL 



(Xa cos 9 + X 4 sin 9) + ^j 2 )+e lKL ( X * 3 sin 9 - xl cos 9) (^ 3 + ^ 4 ) , 

(38) 



{2g+AiK)L _ i 

A = IB 2g + 4tR (Xl cos 9 + xl sin 9) (xl sin 9 - xl cos 9) 2 , (39) 
■02 = 7 — — — (Xi cos ^+X2 sin 9) A \xi sin# - X2 cos#| 2 +|xi cos# + X2 sin 6*| 2 



2g 



(40) 



e (2g-4iK)L _ 1 

•03 = — — — _ — (xi sin 9 - X2 cos^)(xi cos# + X2 sin#) 2 , (41) 



2g - AiK 
e 2gL _ 1 

i>A = 7 Yg S[n9 ~ X l COS ' 



1 2 I 1 2 

A |xi cos 6* + X2 sin sin# — X2 cos9\ 

(42) 
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